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ABSTRACT 

The multifrequency radio continuum and 21cm Hi observations of five blue 
compact dwarf (BCD) galaxies, Mrk 104, Mrk 108, Mrk 1039, Mrk 1069 and I 
Zw 97 using the Giant Meterwave Radio Telescope (GMRT) are presented here. 
Radio continuum emission at 610 MHz and 325 MHz is detected from all the 
observed galaxies whereas only a few are detected at 240 MHz. In our sample, 
three galaxies (Mrk 104, Mrk 108 and Mrk 1039) are members of groups and two 
galaxies (Mrk 1069 and I Zw 97) are isolated galaxies. The radio emission from 
Mrk 104 and Mrk 108 is seen to encompass the entire optical galaxy whereas 
the radio emission from Mrk 1039, Mrk 1069, I Zw 97 is confined to massive 
Hll regions. This, we suggest, indicates that the star formation in the latter 
group of galaxies has recently been triggered and that the environment in which 
the galaxy is evolving plays a role. Star formation rates (SFR) calculated from 
610 MHz emission is in the range 0.01 — 0.1 Mq yr^^; this is similar to the 
SFR obtained for individual star forming regions in BCDs. The integrated radio 
spectra of four galaxies are modelled over the frequency range where data is 
available. We find that two of the galaxies Mrk 1069 and Mrk 1039, show a 
turnover at low frequencies which is well fitted by free-free absorption whereas 
the other two galaxies, Mrk 104 and Mrk 108, show a power law at the lowest 
GMRT frequencies. The fiatter spectrum, localized star formation and radio 
continuum in isolated galaxies lend support to stochastic self-propagating star 
formation (SSPSF). The Hi observations of four galaxies Mrk 104, Mrk 108, Mrk 
1039 and Mrk 1069 show extended disks as large as ~ 1.1 — 6 times the optical 
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size. All the observed BCDs (except Mrk 104) show rotating disk with a half 
power width of ~ 50 — 124 km s~^. Solid body rotation is common in our sample. 
We note that the tidal dwarf (TD) origin is possible for two of the BCDs in our 
sample. 

Subject headings: galaxies: dwarfs, Blue Compact Dwarf (BCD); individual: Mrk 
104, Mrk 108, Mrk 1039, Mrk 1069, I Zw 97 



1. Introduction 



Blue compact dwarf galaxies (BCDs) are star forming dwarf galaxies whose bluer colours 
are attributed to ongoing star formation. They are gas-rich, compact, lo w luminosity 
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■17 to —14) objects with low metal abundances {tqZq < Z < \Zq: 



Izotov et al. 



2006 ). They are not forming stars for the first time, as was predicted earlier ( ISearle et al. 



Searle fc Sargentlll972l ). All BCDs possess a faint low surface brightness component 



that is detected bot h in the optical and in the IR (jCaon et al.ll2005l and references therein; 
Ramya et al.ll2009bl ). implying the presence of old red stars. By analyzing colour maps and 
surface brightness profiles of the low surface brightness (LSB) component, the ages and chem- 
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ical abundances of the underlying host gal axies have been determined (iNoeske et al 
Papaderos et al.l l2002l : iRamya et al.l l2009b| ) . The enrichment of t he ISM in dwarf galaxies 



mainly occurs during these short starburst events (lLegrandll2000[ ). iKunth fc Sargentl (119861 ) 
proposed that if the metals produced during a starburst are immediately mixed with the sur- 
rounding Hll regions, the metallicity will rise very quickly to values of the order of l/50th 
of the solar value which explains why no g alaxy with met allicity lower than I Zw 18 (l/50th 
Zq), a nearby BCD, has ever been found (lLegrandll2000l ). 



A search for quiescent BCDs (QBCDs), carried out by ISanchez Almeida et al.l (|2008[ ) 
indicates that after their bursting phase of a few 10 Myr to a few 100 Myr, BCDs enter the 
quiescent stage. BCDs spend about 30 times more t ime in the quiescent phase. However, 
QBCDs are found to be more metal rich than BCDs (ISanchez Almeida et al.l 120081 ) which is 
yet to be understood. It is notice d that, none of the dwarfs or low surface brightness (LSB) 
galaxies show a SFR equal to zero (ILegrandl (120001 ) and refere nces there i n) wh ich implies even 
during quiescence star formation occurs at a very low rate. ILegrandl (120001 ) has concluded 
through his modelling that the observed oxygen and carbon abundances in I Zw 18 can be 
reproduced by a continuous SFR of 10~^ Mq yr~^after 14 Gyr. However to reproduce the 
present colours, they had to include a bursting episode. All this suggests the existence of 
a weak but continuous regime of star formation in these galaxies. A study of extremely 
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isolated BCDs by IZitrin et al.l (|2009[ ) concludes that the galaxy colours are better explained 
by the combination of a continuous star formation process with a recent instantaneous star 
burst, than by a combination of several instantaneous bursts as suggested previously. 

A few BCDs also emit high ionization lines of He ll A4686 A, [Ne v ] AA3346, 3426 A, 
[Fev ]A4227 A, [Fevil ] AA5146,5177 A, along with broad emission lines of H/3, [Olll ] 
AA4959, 5007 A and Ha in very low metallicity and dense interstellar medium which are 
believed to be due to supernovae (SNe) events and/or stellar winds (llzotov et al.l 120071 . 
and references therein). Whether these low- metallicity BCDs can host an active nucleus 
is another current area of research. Chand ra X-ray observatio ns of star bursting dwarf 
galaxies, such as NGC 1569 and NGC 3077 ( IGrimes et al.ll2005l ) seem to indicate that the 
material is blown out into the halo and consequently even removed from the galaxy leading 
to enrichment of the intergalactic medium. 



BCDs harbour appreciable amounts of dust (IThuan et al.lll999br). confirme d from the 



far-IR (FIR) emission at 60 /zm, 90/im and 140/im ( iHirashita fc Ichikawa 12009). The op- 
tical properties of dust are similar to the dust in the Milky Way ( iHirashita fc Ichikawa 
20091 ). However the dust in BCDs appears to be warmer. Weak polycyclic aromatic hydro- 
carbons (PAH) emission in the bands at 6.2, 7.7, 11.2 and 12.8yum is detected in some 
BCDs. PAH emission is suppressed in most metal-poor B CDs, belieyed to be because 



of a metallicity threshold below wh ich PAHs cease to form (jWu et al.l l2009l : iDwek et al. 



20051 ). lEngelbracht et al.l (120051 ) and lEngelbracht et al.l ( 120081 ) found an anticorrelation be- 
tween the dust temperature and metallicity implying warmer dust at lower metallicities of 
log(0/H) + 12~ 8 and temperature continues to fall with further reductions in metallicity. 
The dependence on metallicity is found out to be ~ down to log(0/H)+12~ 8. The 

change in dust behaviour in terms of PAH emission, FIR colour temperatures and dust /gas 
mass ratio, all near metallicity log(0/H) + 12= 8 indicate that near this metallicity there is a 
general modification of the compo nents of the interstellar dust that dominates the infrared 
emission (lEngelbracht et al.ll2008l ). 



Radio observations which include the 21cm spectral line of H I and radio continuum 
emission are useful in estimating the neutral gas content and kinematics, star formation 
rates and possible signatures of interactions. A sample of BCDs observed in Hi confirrn s 
that metal poor systems tend to be gas-rich low-luminosity galaxies (IHuchtmeier et al.ll2007n. 
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spectrum is attributed t o star forrnation . The FIR-radio correlation of BCDs is similar to 
that of normal galaxies (lYun et al.l 1200 ll ) . 



The initial triggering mechanism, evolution of starburst and evolution of BCDs as a 
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whole is not yet understood. Several mechanisms have been proposed, ranging from internal 
instabilities to external (especially tidal) triggers. If systems are isolated, star formation can 
be explaine d using the stochast ic self propagating star formation mechanism (SSPSF), first 
pro posed by 



Gero. 



ies ( iNoeske et al. 



a et al.l ( 119801 ). Recent studies of large samples of star- forming dwarf galax- 
I2OOII ) that look for faint companions support the hypothesis of interaction- 
induced star formation in BCDs. A lower limit for the fraction of star forming dwarf galaxies 
found with companions is ~ 30% (INoeske et al.ll200ll ). Thus, both the mechanisms are plau- 
sible and it is difficult to quantify the relative influence of the two mechanisms at a particular 
epoch. 

The flve blue compact dwarf galaxies studied here (Mrk 104, Mrk 108, M rk 1039, Mrk 



2009b 




Ramva et al. 


2010) 



In this paper, we present the 21cm Hi tracing the neutral 
atomic gas and radio continuum observations at 240, 325 and 610 MHz tracing the combined 
distribution of thermal and non-thermal radiation for the flve BCDs. This is the flrst time 
that many of these galaxies have been detected at frequencies < 1 GHz. Combined with 
the higher frequency observations from literature, where available, the radio spectra can be 
modelled. The distance to these galaxies is between 20 and 40 Mpc. Table [??] lists the 
general properties of these galaxies. 

Mrk 104 belongs to a loose group UZC-CG 94 consisting of 3 members, with the closest 
member, UGC 4906 an Sa galaxy separated from Mrk 104 in the sky plane by ~ 330 kpc. 
The other member is P GC 26253. Mrk 104 has beeri c lassifled as ha v ing a d ouble nucleus in 
the process of merging (IMazzarella fc Borosoru 119931 ). iRamya et al.l (j2009bl ) resolve the two 
nuclei and note that both show Hll region like spectra thus ruling out the presence of an 
AGN in the centre of the galaxy. Mrk 104 is situated at a distance of 31.2 Mpc and at this 
distance 1" corresponds to ~ 151 pc. 

Mrk 108 classifled as an 10 in NED is one of the four members of the group Holm- 
berg 124. NGC 2820 (of Hubble type SBc) is the closest neighbour. There is also a radio 
bridge connecting NG C 2820 and Mrk 108 to the third member of the group NGC 2814 
( iKantharia et al.l 120051 ) clearly indicating a tidal intera ction. Several other sig natures of hy- 
drodynamic processes are also observed in the group (IKantharia et al.ll2005l ). This galaxy 
hosted a type Hp supernova, SN 1998bm indicating recent massive star production. Mrk 
108 is situated at a distance of 22.2 Mpc and at this distance 1" corresponds to ~ 108 pc. 



Mrk 1039 is a mem ber of the gro up USGC S087 (iRamella et all 120021 ) and LGG 59 
(Lyon group of galaxies: ICarcial 09931 ). DDO 023, DDO 020 and Mrk 1042 are the other 
members of the group. All the companion members are dwarf galax ies. Though locate d 
close to the Eridanus supergroup, it is not considered to be part of it ( iBrough et al.ll2006l ). 
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A type II supernova, SN 1985S is recorded in Mrk 1039. The galaxy is located at a distance 
of 28.8 Mpc (galactocentric distance taken from NED) and at this distance 1" corresponds 
to ~ 139 pc. 

The fourth galaxy in our sample, Mrk 1069 also lies close to the Eridanus supergroup 
but is not considered to be member of the group. No group membership is assigned to Mrk 
1069. This galaxy is located at a distance of 20.7 Mpc from us and at this distance 1" 
corresponds to a distance of ~ 100 pc. 



I Zw 97 is an isolated galaxy with no neighbour within about 50' (525 kpc). iThuan fc Martin 



(Il98ll ) do not detect this galaxy in H I and conclude that the atomic gas surface density is 
< 2.7 X lO^MgMpc-^ and the upper limit on the Hi mass is 7 3 x 10^ M^. Type II SN 



2008bx was discovered recently in this galaxy, and iRamya et al.l (l2009al ) reported the de- 
tection of radio continuum emission at 610 MHz from this supernova. The galaxy is at a 
distance of 36.1 Mpc. At this distance 1" corresponds to ~ 175 pc. 



This paper is structured as follows. Section [2] gives an account of the observations 
and data reduction. Section [3] gives a note on individual galaxies. Section H] is a detailed 
discussion on our results and section [5] summarizes the study. 



Table 1 : : The general parameters of the five galaxies collected from literature. 



Parameter 






Galaxy 








Mrk 104 


Mrk 108 


Mrk 1039 


Mrk 1069 


I Zw 97 


Hubble type'^ 


Pec 


10 pec 


Sc eHtre-on TTtt 


Sa* 




Helio veP (km s~^) 


2235 


1534 


2111 


1562 


2518 


Central vel^ (km s^"'") 


2203 


1574 


2098'' 


1562"^ 


2530 


Group 


UZC-CG 94'' 


Holm 124 


USGC 8087'= 






rnprnbpT's 


UCC 490fi 


NCC 2820 


DDO 023 


UGGA 052 






PGC 26253 


NGC 2814 


DDO 020 












Mrk 1042 






Single dish Hi mass"^ (Mr?,) 


5.0 X 10* 


9.4 X 10^ 


1.2 X 10^ 


8.5 X 10* 


< 7.3 X 10* 


50% line width"^ (km s-^) 


163 


324 


149 


106 




Galactocentric Distance'^-'' (Mpc) 


31.2 


22.2 


28.8 


20.7 


36.1 


linear scale - 1" '^•'^ (pc) 


151 


108 


139 


100 


175 


m.B^ (mag) 


15.1 ±0.4 


15.5 ±0.3 


13.94f 


14.559 


14.9 ±0.3 


Mb (mag) 


-17.46 


-16.47 


-18.46 


-17.32 


-16.76 


Lb (10^ X Lq) 


1.5 


0.6 


3.8 


1.3 


0.79 


Lfir (10^ X Lq) 


0.54 




1.96 


0.88 


0.82 


SN recorded 




SN1998bm 


SN1985S 




SN2008bx 



Thuan fc Martini (11981). - Focardi fc Kelm 
de Vaucouleurs et al.l (|l99ll ). ^ - iDovle et al. 



(I2OO2I) - iRamella et"al] (|2002l ). - Ixhuan et all|l999al) . - NASA Extragalactic 
|2005l ). - Assuming Ho = 73 km s'^ Mpc-\ * - Hyperleda. 
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2. Observations and data reduction 

Hi and radio continuum observations of Mrk 104, Mrk 108, Mrk 1039, Mrk 1069 and 



I Zw 97 were obtained using the Giant Meterwave Radio Telescope (GMRT: ISwarup et al. 



19911 ). GMRT consists of 30 antennas of 45m diameter distributed along a Y. 14 antennas 
are situated in the central 1 km region whereas the rest are spread over a 25 km region. 
GMRT operates at five frequency bands namely 150 MHz, 240 MHz, 325 MHz, 610 MHz 
and 1420 MHz. We observed all the galaxies except Mrk 108 in the dual frequency mode 
in which data in two frequency bands namely 240 MHz and 610 MHz are simultaneously 
observed. We obtained the observations of Mrk 104, Mrk 108, Mrk 1039 and Mrk 1069 at 325 
MHz bands also. All the galaxies except I Zw 97 were observed in the 21cm line of H I in the 
1420 MHz band. We used a bandwidth of 4 MHz for all the galaxies and with 128 spectral 
channels this resulted in a channel width of 31.2 kHz (~ 6.6 km s~^). Appropriate phase and 
amplitude (which also doubled as bandpass) calibrators were selected for the observations. 
Data reduction was carried out using the standard tasks in 

To summarise our data reduction procedure, we imported the FITS file of the raw 
data to AIPS and selected a spectral channel after examining the data on the amplitude 
calibrator which was relatively free of radio frequency interference. This spectral channel of 
the amplitude calibrator was then gain calibrated and bad data were edited. Once the gain 
calibration was satisfactory, the data were used to obtain the bandpass calibration tables. 
In the next step, the visibility data on all sources were bandpass calibrated and then every 
10 channels of the central 100 channels were averaged to avoid the effects of bandwidth 
smearing in the outer parts of the primary beam. The data on both the amplitude and 
phase calibrator of one of the 10 resultant channels were then used to obtain the complex 
antenna gains. After an iterative procedure involving calibrating the calibrator data and 
removing bad data, the final gain tables were generated and the target source data were 
calibrated. 

The data were then imaged by using 9 facets across the primary beam with cellsize of 
1" at 1420 MHz, 25 facets with cellsize of 2" at 610 MHz and 49 facets with cellsize of 4", 3" 
at 240 MHz a nd 325 MHz respectively . The analysed data on Mrk 108 were obtained from 



the authors of iKantharia et al.l (120051 ). The continuum data were then self-calibrated. At 
least three iterations of phase self-calibration were required. After the first self-calibration 
run, it was noticed that the 610, 325 and 240 MHz images improved considerably. A final 
round of amplitude and phase self-calibration was done for each of the data sets. The rms 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under co-operative agreement by Associated Universities, Inc. 
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noise in the final 610 MHz maps is in the range 65-150 /iJy whereas it is 1.2 — 1.9 mJy 
for the 240 MHz images. The rms noise in 325 MHz maps varies from 300-900 fiJj. The 
rms noise in the 1420 MHz continuum maps range from 50 yuJy for VLA Archive image to 
150 yuJy for the maps obtained from GMRT. Self-cahbration at 1.4 GHz did not improve the 
image quahty and hence the iteration procedure was stopped after two rounds of phase self 
calibration. The line-free channels were used to obtain the continuum images at 1.4 GHz. 

The spectral line cube was made after removing the continuum emission using the task 
UVSUB. We also constructed image cubes of angular resolutions varying from around 40" to 
10" to examine both large scale and finer features in the H I distribution. The moment maps 
of these galaxies were created using the momnt task in AlPS. The galaxy I Zw 97 was not 



obser ved in H I as this galaxy was not detected in single dish observations (iThuan Sz Martin 



19811 ). More details of these observations and the results are presented in the Tables |5] and 
|3l We also made images of Mrk 1039 using archival data from the Very Large Array (VLA) 
at 1.4 GHz, 4.8 GHz and 14.9 GHz. 



While all the observed galaxies were detected at 325 MHz, 610 MHz and 1420 MHz, 
only Mrk 104 and Mrk 1069 were detected at 240 MHz. All the observed galaxies were also 
detected in the 21cm spectral line of hydrogen. The contours of the continuum images and 
Hi column density, velocity maps of these galaxies overplotted on the optical images from 
DSS are shown in Figures [iHSl Figure [1] shows the radio Hi line and continuum data maps 
of Mrk 104, Figure [2] displays the radio maps of Mrk 108 at H I line and radio continuum 
bands. Figures |31 ID and E] show the radio continuum and Hi line maps of Mrk 1039, Mrk 
1069 and I Zw 97 respectively. Details of the fiux densities along with rms noise of the maps 
are presented in Table [2l Table [3] shows the Hi properties of these galaxies and Table [5] gives 
the radio-FIR properties of these BCDs. 



3. Results: Notes on individual sources 

1. Mrk 104 : Figured] (a) and (b) show the Hi column density and velocity field of Mrk 
104. The Hi extent (-Djjj, diameter of Hi upto a column density of 5 x 10^^ cm~^) is 
about ~ 2.7 times the optical size (-D25, diameter at 25 mag arcsec"^) of the galaxy. A 
small offset is noticed in the kinematic centre with respect to the optical centre. The 
Hi mass is estimated to be ~ 2.2 x 10^ Mq (refer Table [3]) which is within a factor 
of 2 of the mass estimated from single dish measurements. Figure [H^b) shows the 
kinematics of the H l gas. We note that Mrk 104 shows a distorted velocity field. The 
velocity is redshifted outwards from the centre of the galaxy (see Figure [T|d) reaching 
velocities ~ 2260 km s~^ in the northern parts and velocities of ~ 2250 km s~^ in the 
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southern parts. We detect a distinct cloud near the northern edge of the galaxy at 
velocities ~ 2211 km s~^ (see Figure [TJ)). This cloud has no optical counterpart. 

Figure Wl^c) shows the Hi colunan dens ity contours overlaid on the Ha image of the 
galaxy taken from iRamya et al.l (j2009bl ). The Hi peak is located in between the two 
H II regions as seen in the Ha image. Figure IH^d) shows the position- velocity curve 
drawn along the major axis of the galaxy. Deviations from simple rotation are clearly 
seen in this diagram. As seen in the lower part of the figure, the edges of the galaxy 
show high velocity whereas the central parts show lower velocities. The northern cloud 
seen near declination 53°27'10" shows a very different velocity compared to the edge 
of the galaxy indicating it is not a part of the galaxy. 



Table 2:: The radio continuum flux density and RMS noise on the GMRT images of the five galaxies. 



Galaxy 


RA 


Dec 


Date of obs 


Bands 


Bandwidth 


on-source time 


Flux density 


RMS noise 




hh:7TiV[i:ss 


dd:TfiTfi:ss 




MHz 


MHz 


hrs 


vfijy 


vnjij /beam 


Mrk 104 


09:16:45.5 


+53:26:35 


16/06/08 


240 


6 


5 


7.33 ± 2.6 


1.66 








12/06/10 


325 


16 


2.5 


3.95 + 0.86 


0.22 








16/06/08 


610 


16 


5 


2.40 + 0.3 


0.12 








21/06/08 


1420'' 


1.6 


7.5 


1.66 ± 0.5 


0.18 


Mrk 108 


09:21:30.1 


+64:14:19 


20/08/02 


325 


16 


5 


16.52 + 3.80 


0.77 








06/09/02 


610 


16 


5 


9 03 + 3 6 


0.31 








16/07/02 


1280 


16 


7.5 


1.95 + 0.6 


0.75 


Mrk 1039 


02:27:32.8 


-10:09:56 


16/06/08 


240 


6 


5 


< 5.00+ 1.7 


1.66 








12/06/10 


325 


16 


2.5 


5.8 + 1.1 


0.90 








16/06/08 


610 


16 


5 


7.85 + 0.4 


0.17 


— 


— 


— 


13/07/09 


1420" 


50 


4 


3.00 + 0.2 


0.05 








14/06/00 


4800° 


50 


4 


1.45 + 0.1 


0.05 








16/07/09 


14900" 


50 


4 


1.26 + 0.4 


0.15 


Mrk 1069 


03:08:19.0 


-13:54:11 


09/06/09 


240 


6 


5 


10.33+1.3 


1.99 








12/06/10 


325 


16 


2.5 


12.86 + 2.73 


1.30 








09/06/09 


610 


16 


5 


13.92 + 0.4 


0.06 








31/05/09 


1420'' 


1.6 


7.5 


9.57+1.9 


0.15 


I Zw 97 


14:54:39.0 


+42:01:25 


04/06/09 


240 


6 


5 


< 4.50 + 1.50 


1.50 








04/06/09 


610 


16 


5 


1.14 + 0.9 


0.06 



" Higher frequency data of Mrk 1039 is obtained from the VLA Archives. 

'' This data is obtained after collapsing the line free channels of H i data observed from the GMRT to obtain the continuum map. RMS 
noise here corresponds to the noise in the 50 channels that are collapsed. 
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(a) 



(b) 



(c) 



(d) 

Fig. 1. — : Mrk 104 : (a) Column density, iV(Hi ) contours at (0.5, 1, 2, 3, 4, 5, 6, 7, 7.8) xlO^o cm-^ 
overlaid on the optical B band DSS image, (b) The velocity ma.p{MOMENT 1 map) with velocity contours 
drawn from 2231, 2236, 2241, 2246, 2251, 2256 km s'^ for this galaxy The central contour in black is 2231 
knis^^. (c) The Hi column density cont ours ((0.5, 2, 4, 6, 7.5) xlO^" cm~^) overlaid on the continuum 



subtracted Ha image ([Ramva et al.ll2009bl ). The peak contour shows that the peak Hi emission lies between 
the two Hii regions, (d) Position- velocity curve along the major axis (5° east of north) of the galaxy. Note 
that the major axis coordinates are plotted along the x-axis. The contours are plotted at (2.5, 3, 4, 5)xl.3 
mjy/beam. Note the similar velocities observed in the north and south of the galaxy and the cloud showing 
a distinct identity in the PV diagram. It could be an infalling cloud. The angular resolution of all the maps 
is 21" X 18". 



The galaxy is detected in the radio continuum in all the four bands — 240, 325, 
610 MHz and 1.4 GHz (Figure [T^,f,g,h). The continuum emission encompasses the 
entire optical galaxy and appears to consist of two radio peaks almost coincident with 
the Ha peaks. The global spectral index, a(325,610) = —0.8 (where Sy oc i/") and 
a(610, 1420) = —0.4 imply dominance of non-thermal emission at the lower frequencies 
and probable contribution from thermal processes at the higher frequencies. Figure 
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(i) 



Fig. 1. — : Mrk 104 Continued: (e) contours (1.3 x (-6, -4, 4, 5, 6) mJy/beam) show 240 MHz and 
the grey scale represents the B band optical image. The resolution of the map is 24" x 12". (f ) contours 
(300 X (-6, -4, 4, 6, 8, 10, 11) ^Jy/beam) show 325 MHz overlaid on the grey scale DSS B band optical 
image. Beam size is 24" x 10". (g) contours (140 x (-6, -4, 4, 6, 9, 10) ^Jy/beam) show 610 MHz and the 
grey scale represents the B band optical image. Beam size is 18" x 15". (h) contours (145 x (-5, -4, 4, 5, 
5.5) /ijy/beam) s how 1.4 GHz continu um map and the grey scale represents the Ha narrow band optical 
image taken from lRamva et al.l (|2009b[ ). The 1.4 GHz continuum emission shows emission peaks coincident 
with the bright H ii regions. Beam size is 8" x 7". (i) 610 MHz continuum map (natural weighted) is 
overlaid on the Ha image taken from iRamva et al.l (|2009bf ). The contours levels drawn are 144x(-6, -4, 4, 
6, 8, 10)/.tJy/beam. In this figure and all the figures to follow, North is up and East is towards left. 



[ ITi) sho ws the 610 MHz map overlayed on the Ha images taken from iRamya et al. 



(j2009b ). The radio peak is coincident with the two bright Hll regions seen in Figure 
[U^c) as well. 
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2. Mrk 108 : This galaxy is in a group, namely Holmberg 124 and group interaction is 
clearly noticeable (refer Figure |2]^ a)). Figure [2]J^a) shows a larger region in 325 MHz 
radio continuum in which the group members NGC 2820 (eastern edge), NGC 2814 
(western edge) are also seen. Figure [2](b) shows the H I intensity map. Figure Wic) 
shows the velocity field of Mrk 108 which is distinct from that of NGC 2820. The 
systemic velocity of the galaxy given in NED is 1562 km s~^. We estimate a Hi mass 
of 1.6 X 10® Mq with an error of ~ 30% due to the large H I disk of its neighbour NGC 
2820. 

The galaxy is detected at 325 MHz, 610 MHz and 1280 MHz (see Figure [2] d,e,f). 
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Fig. 2. — : Mrk 108 : Figure (a) shows the triplet in the group Holmberg 124. Contours are 325 MHz 
contours overlaid on B band image. The cross marks the position of Mrk 108. The box drawn in the image 
represents the sizes of the image in Figures (b), (d), (e), (f). Figure (b) represents Hi column density 
contours of Mrk 108 overlaid on optical DSS B band image. The edge on galaxy located on the eastern side 
is NGC 2820 (refer Figure (a)). Contours are plotted at (3, 5, 7, 9, 11) xlO^" cm~^. The angular resolution 
of the map is 23" x 21". Figure (c) shows the MOMENT 1 map including the velocity fields of NGC 2820 
(eastern edge) and NGC 2814 (western edge). Contours are plotted at 1390, 1400, 1410, 1420, 1430 km s^^. 
The kinematical l y dist inct nature of Mrk 108 compared to NGC 2820 is visible. This map is taken from 
Kantharia et al. ( 20051 ) and has an angular resolution of 47" x 34". (d) shows the 325 MHz radio image of 
Mrk 108 zoomed in and overlaid on the optical B band image. Contour levels are 656 x (-4, -3, 3, 4, 6, 8, 10) 
/iJy/beam. The resolution of the map is 22" x 15". (e) shows the 610 MHz contours of Mrk 108 overlaid 
on the optical B band image. Contour levels are 245 x (-6, -4, 4, 6, 8, 10, 12) /^Jy/beam. The resolution of 
the map is 22" x 15". (f) shows the 1280 MHz contours overlaid on optical DSS i mage. Levels are 9 x ( -6, 
-3, 3, 4, 6, 8) /iJy/beam. Resolution of the maps is 6" x 5". Maps are taken from iKantharia et al.l (|2005l ). 



Table 3:: The observation and results of Hi for the four galaxies observed using GMRT. 



parameter 




GalELxy 




Mrk 104 


Mrk 108 Mrk 1039 Mrk 1069 



Date of obs 


21/06/2008 


28/10/2002 


13/07/2009 


31/05/2009 


Bandwidth (MHz/km s j 


4/858 


4/858 


4/858 


4/858 


Cnannelwidtn (KHz/km s j 


0I.25/0.7 


31.25/6.7 


31.25/0.7 


31.25/0.7 


on-source time (hr) 


7.5 


7.5 


7.5 


7.5 


Hi highest resolution (arcsec) 


01// ■. , A Cill 

21 X 18 


11 X 10 


10// , , 1 r\// 

13 X 10 


21 X 13 


RMS in the highest resolution map (mJy) 


1.1 


1.2 


1.9 


2.1 


line width at 50% peak 


68.5 ± 1.6 


51.6 ± 1.9 


124.0 ± 1.1 


61.5 ±0.1 


line width at 20% peak 


152.6 


190.4 


197.5 


104.8 


Total flux (Jy kms^^) 


0.97 ±0.05 


1.42 ±0.41 


6.53 ±0.10 


2.75 ±0.20 


^Hl /^25 ' 


2.66 


2.13 


1.1 


5.43 


M(Hi ) (Afo) 


2.17 X 10^ 


1.6 X 10** 


1.25 X 10^ 


2.71 X 10** 


M(*) (Mq) - 


1.4 X 10^ 


3.1 X 10** 


2.2 X 10^ 


2.13 X lO'' 


M(dyn) {MqY 


1.17 X 10^ 


9.4 X 10** 


4.7 X 10^ 


5.3 X 10^ 


M(Hi )/Lb (Mq/Lq) 


0.15 


0.27 


0.30 


0.20 


M(Hi )/Lk (Mq/Lq) 


0.12 


0.41 


0.40 


0.10 


M(dyn)/LK (Mq/Lq) 


0.83 


2.86 


1.47 


1.99 



" - gives the details of highest resolution map that we have used. 

^ - Except for the galaxy Mrk 1069, the optical diameter {D25) is measured at 25 mag arcsec^^ from lde Vaucouleurs et al.l (|l99lh . For Mrk 
1069, the optical diameter is taken from K band image as gi ven in NED. 
- We adopt the relation M^,/Lk ^ 0.8 (M/Lk)o given bv iKirbv et al.l (|2008l ) for galaxies with B 



i? - 2 ( Vaduvescu et al.l 120071 ) . the 



average colour of the dwarfs. The K band luminosity is es timated from 2 MASS m agnitudes. 

- This is lower by an order of magnitude compared to iThuan fc Martini (|l98lh who used single dish observations and hence included 
contribution from the massive neighbouring spiral NGC 2820. Our interferometric estimate is a better estimate of the mass of this galaxy. 

- M(dyn) is estimated using the simple assumption of circular orbits i.e. using M — 
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3. Mrk 1039 : The Hi intensity and kinematics are shown in Figure |3]J^a),(b). A warped 
gas disk is visible in Figure whereas the stellar disk does not show any signature 
of the warp. The velocity field (Figure [3I^b)) shows the rotating gas disk with small 
distortions in the out er parts espec i ally in the western half. Such small kinks and bends 
were also noticed by iThuan et al.l ( l2004j ). The velocity dispersion of Mrk 1039 varies 
from ~ 15 km s~^ at the edges of the galaxy to 45 km s~^ in the central regions. The 
Hi mass is estimated to be 1.25 x 10^ Mq (Table [3]) which is in good agreement with 
the mass estimated from single dish measurements. The Hi disk is about ~ 1.1 times 
the (-D25) optical size of the galaxy. Figure [3](c) shows the Hi column density map 
at a resolution of 13" x 10". The smooth H I distribution seen in the lower resolution 
image is resolved into several high column density clumps, noticeably the clouds on 
the eastern and western edges. Figure El^d) shows the position- velocity diagram along 
the major axis (90° east of north). The disk seems to execute solid body rotation till 
the outskirts of the galaxy. Figure [3]^e) shows the H I column density (high resolution) 
contours overlaid on the Ha images taken from iRamya et al.l (120101 ). It is seen that 
the high column density Hi gas overlaps with the strong Hll regions in the eastern 
side of the galaxy, from where radio continuum emission is also detected. 

The continuum maps at 325, 610 MHz (GMRT data), 1.4 and 14.9 GHz (from VIA 
archival data) are shown in the Figures [Hf-i. The radio emission appears to be coupled 
to the bright Hll region near the eastern edge of the galaxy. Note that the emission 
at 610 MHz is more extended compared to that at 14.9 GHz indicating the presence of 
non-thermal emission. Fi gures [3](i) and (k) sh ow 325 and 610 MHz contours overlaid 
on the Ha image shown in lRamya et al.l (120101 ). This galaxy hosted a type H supernova 
SN 1985S close to the region from where we detect radio continuum emission (marked 
by a star in Figure [3]). We do not detect the galaxy at 240 MHz with a 3a limit of 5 
mJy. The integrated spectral index between 610 MHz and 1.4 GHz is a ~ —0.8 and 
a(325,610) = +0.48. 

4. Mrk 1069 : The Hi momen t maps are shown in Figure |H^a),(b). This galaxy, classified 
Sa type by Hyper ledaB (jPaturel et al. 2003) shows a large H I disk (see Figure Hl^a)) 



as a I 



about 6 times (Table [3]) the optical size of the galaxy. The H I disk shows a slight warp 
in the northern parts. The southern part of the disk appears to be abruptly truncated. 
The galaxy shows rotation (Figure Hl^b)). A slight offset between the optical centre 
and the kinematic centre is noticed in this galaxy similar to Mrk 104. We estimate a 
H I mass of 2.7 x 10^ Mq (Table [3]) which is lower by a factor of about 3 compared to the 
single dish estimate. Higher resolution (21" x 13") column density distribution is shown 



^http:/ /leda. univ-lyonl.fr 
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Fig. 3. — : Mrk 1039 : (a) The Hi column density contours at (2, 4, 6, 8, 10, 12, 14, 16) xlO^" cni'^ 
are overlaid on DSS B band image. Angular resolution of the map is 48" x 31". The Hi disk is warped. 
The Hi emission (-Djj j) is ~ 1.1 times the optical size {p^^ of the galaxy. The box drawn in the picture 
represents the size of the images in Figures (c)-(k). (b) The velocity field of Mrk 1039 is shown in this figure. 
Contours are plotted at 2020, 2040, 2060, 2080, 2100, 2120 km s'^ (c) The high resolution (13" x 10") 
Hi column density map of Mrk 1039. The Hi gas resolves into several clouds at the higher resolution. 
Contours are drawn at 4, 8, 12, 16, 20 xl0^° cm~^. (d) Position-velocity curve along the major axis (90° 
east of north) of the galaxy. The contours are plotted at 3.1x (4, 5.5, 6.5, 8, 9) mJy/beam. This P-V diagram 
is created using the lowest resolution (48" x 31") Hi cube and shows the solid body rotation in the galaxy. 



(e) T he high resolution Hi column density contours are overlaid on the Ha image taken from lRamva et al 



( 20101 ). The levels are plotted at (5, 12, 20) xlO^ 



The highest column density (2 xlO^^ cm ^) Hi is 



coincident with the strong Ha emitting Hii region in the east. 



in Figure 111(c). Figu re HTd) shows t he hi gher resolution column density map overlaid 
on the Ha images of iRamya et al.l (120101 ). There appears to be an offset between the 
higher column density H I peaks and the Ha peaks. The position- velocity curve (refer 
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Fig. 3. — : Mrk 1039 Continued: Figure (f) shows 325 MHz contours overlaid on the optical DSS image. 
Contour levels are 1.2 x (-4, -3, 3, 4, 6, 8) mJy/beam. The resolution of the map is 25" x 17". (g) shows 
610 MHz contours overlaid on the optical DSS image. Contour levels are 165 x (-6, -4, 4, 6, 8, 10, 12, 14, 
18) /iJy/beam. The resolution of the map is 8" x 6". (h) shows the 1.4 GHz contours obtained from VLA 
archival data overlaid on the optical DSS image. Levels are 84.0 x (-8, -4, 4, 8, 12, 16, 20, 24) /xJy/beam. 
The resolution of the map is 8" x 7". (i) shows the 14.9 GHz contours obtained from VLA archival data 
overlaid on the optical DSS image. Levels are 207.0 x (-4, -3, 3, 4, 5) /zJy/beam. The emission is localised 
being coupled to a star- forming region in the east. Note that a Hi cloud is also seen to be coincident with 
this region. Angular reso l ution of this map is 8" x 7". (j) Shows the 610 MHz contours overlaid on Ha image 



taken from 



Ramva et al. ( 2010l) . The contour levels are 1 65 x (-6, -4, 4 6, 8, 10, 12) /xJy/beam. (k) Shows 



the 325 MHz contours overlaid on Ha image taken from Ramva et al. ( 2010l ). The contour levels are 1.2 x 
(-4, -3, 3, 5, 6, 8) mJy/beam. The radio emission encompasses all the Hii regions seen in the Ha image. 
The star marked in all the images towards south-west direction reperesents the location of the SN 1985S. 



Figure m^e)) along the major axis of the galaxy shows rotation. 

Radio continuum emission from Mrk 1069 is detected at 325 MHz, 610 MHz and 1420 
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Fig. 4. — : Mrk 1069: Figure (a) shows Hi column density contours (1, 2, 3, 4, 5, 6, 7, 8x 10^" cm"^) 
overlaid on the B band DSS image. Angular resolution of the map is 41" x 35". Note the truncation of the 
H I disk in the southern part of the galaxy and the large H i disk compared to the optical disk. The box 
represents the size of the images in Figures (c)-(j). (b) Hi iso-velocity map [MOMNT 1). Contours are 
plotted between 1515 and 1615 km s"^ in steps of 10 km s~^. The velocity field shows a disturbed rotation 
field, (c) Higher resolution (21" x 13") Hi column density map for this galaxy is shown here. Contours are 
plotted at 6, 8, 10 xlO^" cm^^. Only the central 1'.3 x 1'.3 of the galaxy shown in (a) is zoomed-in and 
shown in this map. Note the fragmentation of the Hi disk into s maller clouds, (d) T he higher resolution 
H I column density map is overlaid on the Ha image taken from iRamva et al.l (j2010l ) . Contour levels are 
same as in (c). It is noticed that the highest column density contour is off-centered from the Ha image, (e) 
Position-velocity curve along the major axis (60° east of north) of the galaxy. The contours are plotted at 
2.3 X (3, 4, 5, 6, 7) mJy/beam. This P-V diagram is created using the lowest resolution (41" x 35") Hi cube. 



MHz (Figure Hf-h respectively). The radio continuum emission arises in the central 
parts of the galaxy and extends to the north. The global spectral index between 610 
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Fig. 4. — : Mrk 1069 Continued : Figure (f) shows 325 MHz contours (levels = 1.7 x (-5, -4, 4, 5, 6, 7, 
8, 9) mJy/beam) overlaid on DSS B band optical image. The resolution of the map is 25" x 16". (g) shows 
610 MHz contours (levels = 70.0 x (-6, -4, 4, 5, 6, 8, 12, 16, 20, 28) ^Jy/beam) overlaid on DSS B band 
optical image. The resolution of the map is 5" x 5". (h) shows the 1.4 GHz continmun map obtained from 
GMRT overlaid on the DSS optical image. Contour levels are 140.0 x (-6, -4, 4, 6, 8, 10) //Jy/beam. The 
resolution of the map is 5" x 5". Note the bipolar feature extending along the minor axis of the galaxy, (i) 
shows the spectral index map in gray scale plotted with 1.4 GHz contours at (-6, -4, 4, 6, 8) x 140/iJy/beam. 
Spectral index varies from —0.8 to +1.5. Angul ar resolution o f the maps is 5" x 5". (j) shows the 610 MHz 
contours overlaid on the Ha image taken from Ramva et al. ( 2010l ). The contour levels are 70.0 x (-6, -4, 
4, 8, 12, 16, 28) /LtJy/beam). The resolution of the radio map is 5" x 5". The radio continuum emission is 
coincident with the Ha emission and has a similar extent. 



MHz and 1.4 GHz is ~ -0.4 and between 325 and 610 MHz is ~ +0.13. 

The 610 MHz and the 1.4 GHz images show a bipolar outflow-hke feature extending 
along the minor axis. Interestingly, this feature shows a slightly flatter spectrum as 
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compared to the surrounding emission as seen in the spectral index map between 610 
and 1420 MHz (Figure lUfi)). W e note that a similar feature has been seen in the galaxy 
F08208+281 6 Jvin et al.lboosh Figure shows 610 MHz contours overlaid on the 



Ha image of iRamya et al.l ( 120101 ). it is noticed that the extent of emission in 610 MHz 
and in the Ha are the same. However, no obvious bi-polar like feature is visible in 
the Ha map. Higher sensitivity and higher resolution data are required to confirm the 
existence and nature of this radio feature. 

An extended source situated at 03^08"^21".7 -13°54'42".l (refer Figure SJ^g)) east 
of Mrk 1069 is detected at 610 MHz with a flux density of 1.43 ± 0.18 mJy. This 
background source is identified with 2MASX J03082177-1354425. 

I Zw 97 : We do not detect the galaxy at 240 MHz upto a 5 a limit of 7.5 mJy on 
the flux density of the galaxy. We detect radio emission from this galaxy at 610 MHz 
(see Figure [5]J^a)) with an integrated flux density of about 1.14 mJy. This is the first 
detection of this galaxy in radio band s to the best of our knowledge. We note that 



it ha s not been detected in the NVSS (jCondon et al.l 119981 ) and FIRST (IBecker et al. 



19951 ) surveys at 1.4 GHz. The emission at 610 MHz arises in two discrete regions of 



the galaxy. We also note that a type II sup ernova SN 2008bx was discovered in this 
galaxy on 2008 April 22 JPuckett et al.ll2008h . 
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Fig. 5. — : I Zw 97 : (a) 610 MHz contours superposed on the DSS B band optical image. Contours 
are plotted at 65.0 x (-3, 3, 4, 6) uJy/beam. The reso lution of the map is 6" x 5". (b) 610 MHz contours 
overlaid on Ha image taken from Ramva et al. ( 2009bl ). The star mar ks the position of th e supernova SN 
2008bx. This galax y is not detected at 1.4 GHz either in the NVSS |Condon et allllQQSl) or the FIRST 
iBecker et al.lll995h surveys. 



The southern component seen in our 610 MHz image is emission associated with the 
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supernova SN 2008bx (iRamya et al.ll2009al ). The radio continuum emis sion coincides 



fairly w ell with the Ha emission (see Figure [5t^b), Ha image taken f rom iRamva et al. 



(2009b|)). The galaxy is not detected in single dish Hi observations (iThuan &: Martin 
198lh . 



Discussion 



4.1. Radio continuum spectra 



Figure |6] shows the observed radio spectra of the five galaxies. The spectra of Mrk 1069 
and Mrk 1039 are seen to turn over at frequencies below 610 MHz, while, the spectra of 
Mrk 104 and Mrk 108 show a power law upto the lowest GMRT frequencies. The observed 
radio continuum emiss ion for Mrk 1039 a nd Mrk 1069 is mostly confined to the star forming 
regions traced by Ha. iKlein et al.l ( 1l99ll ) noted that the break frequency is below 1 GHz in 
their sample of BCDs. We note that Mrk 108 (dotted line in Figure [6]) is projected close 
to NGC 2820 and hence it is difficult to separate the emission of NGC 2820 from Mrk 108 
at the low frequencies where the beamsize is larger. There is a possibility that the lower 
frequency emission (i.e. frequencies less than 1 GHz) includes contribution from the spiral 
companion. 

We have combined VLA archival data at higher frequencies with our low frequency 
GMRT data for Mrk 1039. The spectrum (Figure El flattens for z/ > 1.4 GHz and the 
emission appears to be dominated by thermal emission. Since the higher frequency data is 
not available for o ther galaxie s in ou r sample, we cannot comment on their behaviour at 
these frequencies. iKlein et al.l ( 1l99ll ) have noted fiat spectrum at higher frequencies with 
non-thermal spectral index ranging from —0.7 to —2 for some g a laxies in their sample. 
However the low frequency observations (325 MHz) by iDeeg et al.l (119931 ) suggest that the 
average spectral index of the nonthermal emissi on is ~ —0 .7. Moreover, the thermal fraction 



at 1 GHz is seen to vary from about 5% to 72% (IKlein et al.lll99l[ ). Our sample i s a subset of 



a larg er opt ically select e d sara ple of BCDs and hence distinct from the sample of iKlein et al. 



a larg 
(1991 


er opx 
) and 
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(1993) 



continuum emission using four models: (1) the observed emission is non-thermal synchrotron 
emission; (2) the observed emission is a combination of thermal free-free and non-thermal 
synchtrotron emission along the line of sight; (3) the non-thermal emission at low frequencies 
is free-free absorbed by the thermal region in the same volume; and (4) the observed emission 
is a combination of both thermal and non-thermal, and the non-thermal emission is free-free 
absorbed by the thermal gas of optical depth ti. 
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freq (MHz) 



Fig. 6. — : The observed spectra of the five galaxies studied here. The points are the observed data and 
we have connected these by a hne to help the reader distinguish the spectra of different galaxies. Solid line 
represents the data for Mrk 104, dotted line Mrk 108, short dashed line Mrk 1039, long dashed line Mrk 
1069 and the dot dashed line I Zw 97. The upper hmits at 240 MHz for Mrk 1039 and I Zw 97 are shown 
as arrows. In the plot, I Zw 97 is shifted by —125 MHz for clear viewing. Note the differences in shapes of 
spectra of different galaxies. 
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Since these cases well represent what we understand about radio emission from galaxies 
and they explain the observed low frequency spectrum of the BCD galaxies, we have not 
explored the case of synchrot r on sel f-absorption which we believe will not be important 



in these galaxies. iDeeg et al.l (Il993[ ) have shown that both free-free absorption at lower 
radio frequencies and time-dependent electron injection models explain the flattening of the 
spectrum at lower radio frequencies. However since we detect a clear low frequency turnover 
in only two galaxies, we have not tried to distinguish between the two models. Figure [7] 
shows the best fits using one of the four models given below to the observed points. 

The four equations which were fitted were: 

(1) The observed emission is non-thermal synchrotron emission which is fitted with 
power law index a: 

Siv) = c,iy" (1) 

(2) The observed emission is a combination of non-thermal (Snth) and thermal free-free 
emission {Sth)'- 

^(z/) = ^„i,(ocz/") + 5i,(ocz/-°-i) (2) 

(3) The low frequency emission is free-free absorbed by the thern ial material i nterrn ixed 



with the emitting region with optical depth ti (Equation taken from iDeeg et al.l (119931 )): 

S{u) = ciZ/(°+'-^) [1 - exp(-riz/-2.i)] ^3) 

(4) The observed emission is a combination of both non-thermal and thermal emission. 
The non-thermal component at lower frequencies is free-free absorbed by the t hermal gas 



i nterrn ixed with the emitting region with optical depth ri (Equation taken from iDeeg et al. 



(Il993l ) and modified to include thermal emission): 

S{u) = CiZ/("+2.i) _ exp(-riZ/-2-i)] + z/-^-^ (4) 

In the above equations, Ci is a constant, while C2 is the thermal emission at frequency 
u. While a power law i.e. model (1) is the best fit to the spectrum of Mrk 108, the observed 
spectrum of Mrk 104 is best fit by including contributions from thermal and non-thermal 
emission i.e. model (2). The observed spectrum of Mrk 1069 was best fit by model (3) i.e. 
free-free absorption of synchrotron emission giving rise to a low frequency turnover. Mrk 
1039 is best fit using model (4). These fits are shown in Figure [7] and the parameters are 
listed in Table HI Since the number of input points are small, the observed points were 
first spline interpolated. The interpolated points were then given as in put data points to the 
programme which used the Levenberg-Marquardt algorithm as given by lPress et al. Jl993h to 
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find the best fit. The fitted parameters and the physical quantities derived using these fitted 
parameters are given in Table |H The reduced values are given in Table HI Galaxy-wise 
results from the modelling procedure are given below: 

Table 4:: Parameters of the model fits to the observed data. 



Parameters Mrk 104 Mrk 108 Mrk 1039 Mrk 1069 

Model eqn 3 eqn 2 eqn 4 eqn f 

X^^grf 0.39 4.91 2.34 1.124 

ci,C2 — — 10.06 ±0.17, 1.41 ±0.03 98.18 ±2.78 

a. -1.67 ±0.26 -1.36 ±0.12 -1.85 ±0.10 -0.51 ±0.12 

n — — 0.42 ±0.02 0.13 ±0.01 

EM from Ti (pc _ _ (1.26 ± 0.10) x 10^ (3.67 ± 0.20) x 10^ 

n(Hii ) cm-3 — — - 354 " - 192 ° 

fraction of '^thiAGH., (%) 80% — 45% — 



° - Assuming a size of 10 pc for the emitting region. 



Mrk 104: The best fit model (Figure [7]) for this galaxy is when the fitted emission 
consists of both thermal and non-thermal emission. This is also obvious from the observed 
spectrum which shows that the emission at 1.4 GHz is more than expected from a power 
law. We estimate the thermal fraction at 1.4 GHz of ~ 80% and the spectral index, a of 
the synchrotron spectrum is —1.67 (Table H]). Compare this with a typical thermal fraction 
of < 10% and synchrotron spectr al index of ~ —0.7 or —0.8 for normal spiral galaxies 
( jPale fc Helod I2OO2I : ICondoru 119921 ) . However, no signature of free- free absorption towards 
the lower frequencies is seen down to 240 MHz. The fraction of thermal emission in the 
galaxy is fairly large implying a young burst of star formation. iRamya et al.l (j2009bl ) show 
that this galaxy has an abundant 4 Gyr population + ~ 500 Myr population and young 
burst of star formation which is ~ 10 — 15 Myr old. Moreover they show that the major 
contribution to the total mass in this galaxy is the Gyr old population whereas the dominant 
contribution to the total light comes from the ~ 500 Myr population. 

Star formation rate (SFR) is calculated using the 610 MHz emission as given by equa- 
tion 21 of ICondoru ( 1l992h . The spectral index a obtained from the fit is used to calculate the 
SFR. The SFR is 0.023 ± 0. 011 Mcr, vr'^for M rk 104 (Table ^. The SFR estimated using 
the FIR luminosity, following Iyuu et al.l j200lh is 0.12 ±0.04 Mq yr^^ The SFR (> 10 Mq) 
estimated from the Ha fiux is 0. 04 M(^^ yr~^ and is in the range 0.0003 — 0.004 Mq yr~^ for 
individual star forming regions ( iRamya et al.ll2009bl ). The SFR calculated using Ha fiux is 
in reas onably good agreeme nt with the SFR found using the 610 MHz emission. The param- 
eter q ( ICondon et al.lll99ll ) which is a measure of the starburst versus AGN nature of the 
radio emission is also estimated and tabulated. Normal galaxies which show the FIR-radio 



- 26 - 



0.1 



0.1 



I I I I I 1 1 — r 

Mrk104 




1 

freq (GHz) 



1 


1 1 1 1 1 M| 1 1 1 1 1 1 M 1 1 


1 1 


- ( 

- ' 


Mrk1039 










1 


1 1 'i 1 1 


1 i' 



1 10 
freq (GHz) 




J I I L. 







4^ 10 

X 
3 



1 III 


1 1 1 1 
MrklOS 


1 III 


1 , , 1 



freq (GHz) 
1 — 



T 1 1 — TT 



Mrk1069 




J I I I I I I 



1 

freq (GHz) 



0.5 
freq (GHz) 

Fig. 7. — : The model fits to the observed spectra of four galaxies. In each of the panels, the dashed line 
represents the thermal component, dotted line represents the non-thermal component (modified by free-free 
absorption in the case of Mrk 1039 and Mrk 1069). The solid line represents the best fit to the observed 
points. The best fit model to the galaxy Mrk 104 is the combination of synchrotron emission and thermal 
bremstrahlung emission. Mrk 108 is best fit with a power law spectrum. Mrk 1039 is best fit assuming a 
combination of thermal and non-thermal emission, and the non-thermal emission is absorbed by the thermal 
gas mixed in the region. Mrk 1069 is best fit with synchrotron emission and thermal free-free absorption at 
lower frequencies. Table |3] gives the values of the fitted parameters. 



(1.4 GHz) correlation have a value of g = 2.35 ± 0.2 flCondonlll992l : lYun et al.l 120011 ) . while 
q is 2.44 ±0.02 (Table [5]) for Mrk 104 which is within the observed scatter for the correlation. 



Mrk 108: The galaxy located in the group Holmberg 124 is best fit (Figure [7]) with a 
power law spectrum. The spectral index of this non-thermal emission is ~ —1.36 (Table H]) 
and we note that the possibility that the lower frequency flux density is overestimated due to 
contribution by the large companion NGC 2820 cannot be ruled out. If the error introduced 
by NGC 2820 is small and the steep index is real then it could be due to galaxy-wide star 
formation induced by tidal interaction. The SFR (refer Table [5]) calculated using the flux 
at 610 MHz is 0.051 ± 0.036 respectively, q value could not be estimated due to the lack of 
IRAS FIR fluxes in the literature. 



Table 5:: Radio-FIR properties for the five BCDs. Icr errors are quoted in the brackets. 
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° - FIR data is o btained from IRAS Faint Source Catalog (source: NED) and the parameters are calculated using the equations given in 
Yun et all (l200ll) . 

^ - The non-th ermal fraction of the total 610 MHz fiux obtained from the model is used to estimate the SFR using the equations given in 
CondonI (Il992l) . 

- The Ha SFR (> 10 Mq) is c alculated using the e quation of Kennicutt ( 1983 ) and quoted from Ramva et al. ( 2009bl) . 

- The Ha SFR is quoted from lRamva et all |2010t ). 
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Mrk 1039: The observed emission of Mrk 1039 for which we have also analysed and 
included higher frequency data from the VLA archives shows a flattening at higher frequen- 
cies and turnover at lower observed frequencies. The observed emission is well fitted by a 
combination of thermal and non-thermal emission and the non-thermal emission is absorbed 
at lower frequencies by thermal gas as shown in Figure [71 The thermal and non-thermal 
emission are mixed. This is possible in a young (~ 10^~^ yr) starburst region which is also 
recording supernovae as the massive stars evolve but the relativistic electrons have not yet 
had time to diffuse through the galaxy. The spectrum is modelled using equation (4). The 
optical depth is ri=0.42. This gives an emission measure (EM) of 1.26 x 10^ pc cm~^ for an 
assumed temperature of 10000 K. The Hll regions electron density is then rig ~ 354 cm~^ 
assuming a 10 pc size for the region. The parameters of the fit are given in Table ID Using 
this and following iHunt et al.l (120041 ) . the total number of ionizing photons is calculated to 
be 9.82 X 10^^ photons/sec which corresponds to 9820 07 V stars in that region. The H/9 
flux calculated using the thermal radi o flux is 4.72 x 10~^^ erg cm~^s~^, which agrees with 



the values quoted by lYin et al.l (120031 ) and iHuang et al.l (119991). This i mplies that there is 
a young and intense star formation in this galaxy, iHuang et al.l (119991 ) have estimated an 
age of 4 Myr for starforming regions in the galaxy. The spectral index of the non-thermal 
emission is fairly steep, a ~ —1.85 (Table Hj). We note that a type II supernova SN1985S 
has occurred to the south west of the centre of the galaxy (shown by a star in Figure |3]). 
Although its location is close to the part of the galaxy from which radio emission is detected, 
no radio emission is seen to be coincident with it. The 610 and 325 MHz emission is clearly 
more extended than the image at 15 GHz obtained using VLA archival data (see Figure 
in]i,j,k) but does not encompass the supernova. 

We estimate the q paramete r to be 2.90 ±0.03 (Table [5]). A similar value for this galaxy 
was also obtained by I Yin et al.l ( 120031 ) (g = 2.94). SFR calculated using FIR luminosity 
is 0.48 ± 0.09 Mq yr-^ Using the 610 MHz emission, we estimate a SFR of 0.060 ± 0.016 
Mq yr-i for this galaxy. The SFR(> 10 Mq) quoted in the table is 0.173 M© yr'^ but the 
SFR for individual knots varies from 0.01-0.15 Mq yr~^. 

Mrk 1069: The observed spectrum of Mrk 1069 shows a clear turnover at lower fre- 
quencies and hence was modelled using a combination of non-thermal emission and free-free 
absorption at the lower frequency. Higher frequency archival data was not available for this 
galaxy. The best fit is shown in Figure |71 The spectral index of the non-thermal emission is 
~ —0.51 and the optical depth of the absorbing thermal gas is 0.13. Assuming an electron 
temperature of 10^ K implies a modest EM of 3.7 x 10^ cm^^ pc (Table H]). For an H ll region 
of size 10 pc, this implies an electron density of about 192 cm~^. The radio emission is in- 
tense near the optical centre of the galaxy and extends towards the north along the major 
axis as seen in our 610 MHz image. We also detect a bipolar outflow like feature extending 
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along the minor axis of the galaxy in our 1.4 GHz image (Figure |l]h) and hints of the same 
are visible at 610 MHz. Similar feature is not detected in Ha (Figures Hlj) but our images are 
not deep enough. The feature needs t o be confirna e d with higher frequency, higher resolution 
radio continuum data. We note that Hunt et al. J2005h have detected a similar outflow in 
the galaxy I Zw 18. We estimate q = 2.32 ± 0.06 (Table [5]) which agrees with the FIR-radio 
correlation within errors. The SFR calculated for Mrk 1069 using the FIR luminosity is 
around 0.19 ± 0.04 Mq yi'^ and from 610 MHz flux 0.103 ± 0.026 M© yr'^ 

I Zw 97: The galaxy, I Zw 97 is detected only at 610 MHz and is undetected in any 
other frequency. Due to this, the modelling of the spectra could not be done. The SFR 
calculated using the FIR luminosity is 0.17 ±0.05 M p) yr~^(Table[5 l ) and u sing 610 MHz flux 
0.023 ± 0.023 Mq jt-\ Using the Ha flux quoted in lRamva et"aLni2009bh . we find the SFR 
is 0.001 — 0.03 Mq yr^^ if we consider the production rate of stars of masses greater than 
1OM0. 

For all the five galaxies, the SFR calculated using FIR emission is higher than SFR 
calculated using the radio continuum emission at 610 MHz by a factor ranging from 2 to 10. 
The SFR calculated using radio emission is comparable to the SFR calculated for individual 
star forming regions from Ha for Mrk 1039 and I Zw 97, whereas, the SFR for Mrk 104 
matches well with the SFR calculated from the total Ha emission. 



We also e xamined th e natu re of the radio spectra of the BCDs s tudie d in the radio 
continuum by iKlein et al.l ( Il99ll ) and iDeeg et al.l ( 119931 ). Klein et al.l ( 1l99ll ) have studied 
23 galaxies an d a smaller sampl e of seven galaxies from this has been studied at lower 
frequencies by iDeeg et al.l ( 119931 ). In our sample of five, we notice that the galaxies in 
groups (i.e. Mrk 104 and Mrk 108) appear to show a power law spectrum to the lowest 
observed frequency whereas the relatively isolated galaxies (Mrk 1069, Mrk 1039) show a 
low frequency turnover. Although Mrk 1039 is in a group, we note that the other members 
are low mass dwarf galaxies unlike in the case of Mrk 104 and M rk 108 whic h have massive 
spirals as members. We examine the behaviour exhibited by the iDeeg et al.l (119931 ) sample. 
We have tabulated the results in Table |6l The first column gives the name of the galaxy, the 
second the group name or notes its isolated status, the third the number of group members. 
The fourth and fifth columns list the radio spectral indices of the galaxy for two sets of 
frequencies, and the last column gives the reference to the information collected. For the 
galaxies whose membership is not determined, its status is kept as 'single ?' and also the 
information regarding the nearby galaxy and possible partner is mentioned, as obtained from 
Hyperleda. Interestingly, the spectra of galaxies which are relatively isola ted show fia t teniri g 
at the lower frequencies (less than 1 GHz) both in our sample and that of iDeeg et al.l (119931 ). 
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Thus, three out of seven galaxies from iDeeg et aLl (119931 ) are isolated, and show flat- 
tening at the lowest frequency (325 MHz) that they observed. The galaxy Mrk 1069 which 
shows flattening at the lower frequency is isolated, while Mrk 1039, also showing flattening, 
is a member of a group having only dwarfs as group members. The flattening is probably 
related to the compact radio emitting, star forming regions. The localised emission could be 
a result of the trigger of star formation being relatively local in the case of is olated BCDs, 



follow ing the scenario of stochastic self propagating star formation (SSPSF; iGerola et al. 



(jl980l )). Thus, from Table Owe can derive two inferences; (a) the overall spectrum of rel- 
atively isolated galaxies is flatter compared to BCDs in denser environments and (b) the 
spectrum turns over at low frequencies (z/ ~ 610 MHz) for relatively isolated galaxies. How- 
ever, these are preliminary results based on a small sample and should be confirmed with a 
study of a larger sample of BCDs in radio. 



Table 6:: Spectral index of a sample of BCDs including our data and Deeg et al. ( 19931) ; D93 data. Columns 
2 and 3 are obtained from NED and Hyperleda. 
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4.2. Atomic gas distribution 



The atomic gas is much more extended than t he stellar extent, a feature generally 
seen in most late type galaxies and also most BCDs (IThuan et al.ll2004j ). The typical peak 
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H I column density in these galaxies is ~ 10^^ cm ^ at a resolution of 45" x 35" ( ~ 5.0 x 5.0 
kpc). Hence most of the galaxy has lower H I column densities. According to ISkillman 



(119871 ). iTaylor et al.l ( 119941 ). the star formation in dwarf galaxies occurs only above a con- 
stant threshold Hi column density of A^(H l) ~ 10^^ cm~^. However, column densities of 
10^*^ c m~^ are fairly corn mon in metal deficient dwarf galaxies ( lEkta et al.ll2006l : lEkta et al. 



20081 : iBegum et al.l |2008| . and refe rences therein). A preliminary study of a small sample of 



galaxies from the FIGGS survey (IBegum et al.ll2006l ) suggested that un like brighter dwarfs , 
the faintest dwarf galaxies do not show well-defined threshold density (IBegum et al.ll2008l ). 
Thus, the requirement for the threshold column density may not be uniform. 

The H I m asses (Table B l) calcul ated here are in the ra n ge 10 '^ — 10^ Mq similar to the 
ones quoted in iThuan et al.l (Il999al ) and iThuan fc Martini (jl98ll ) . A large discrepancy in 
Mrk 108 is noticed. Since Mrk 108 is located on the edges of NGC 2820, our mass estimate 
is more accurate. 

Figure [8] shows the observed integrated line profiles of H I emission for all the four 
galaxies. Wide double-horned profiles typical of edge-on rotating disks are observed for Mrk 
1039 and Mrk 1069. The widest line is seen from Mrk 1039 with a width at 50% of peak 
of about 124 km s~^. We fit the observed profiles with gaussian components. The Hi line 
profile (Figure [8]) of Mrk 104 was best fit with two components centered at velocities 2211 
and 2279 km s~^ with widths 57 and 69 km s~^, respectively. The widths estimated here 
are smaller than the w idths quoted in si n gle dis h observations (see Tables [??] and [3]) by 
Thuan fc Martini (1198 ll ) and iThuan et al.l (Il999al ). The component at 2211 km s~^ is from 
the Hi cloud seen to the north of Mrk 104 which has no optical counterpart. The velocity 
field of Mrk 104 is distorted as shown in Figure [1] (b) — we also note that it is classified by 
NED as a peculiar galaxy. The velocity field is receding on either side of the centre. The 
Hi distribution across the galaxy shows a steep rise in the column density on the eastern 
side of the galaxy and diffused gas in the SW and NW. Note that the closest projected part 
of the galaxy to this cloud has velocities of ~ 2260 km s~^. Thus either (1) the cloud is 
falling onto the galaxy from behind; a remnant of tidal interaction in the group or (2) the 
cloud has been pulled out of the galaxy. 

The nearest neig hbour to Mrk 104 is the galaxy UGC 4906 whose integrated H I profile 
(ISpringob et~aDl2005h shows it to be a rotating mildly lop-sided Hi disk. PGC 26253 (2297 
kms~^) ~ 17'. 5 (159 kpc) away towards south of Mrk 104 is the other member in the 
group. A separate cloud-like entity is detected to the north of Mrk 104. The galaxy has a 
distorted velocity field. All of the above seem to suggest that Mrk 104 has experienced tidal 
forces. Using the peak velocity from the P-V diagram, we estimated the dynamical mass 
using iy'^R/G). The M(Hl )/Lk and M(dyn)/Li^ ratios calculated are given in Tabled 
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M(Hl )/Lk is about Q.12Mq/Lq and M{dy\i)/LK is about O.83M0/L0 for this galaxy and 
the width of the Hi hne is around 69 km s~^. Most of the parameters seem to suggest that 
the dark matter content is low and a tidal dwarf origin for Mrk 104 is possible. 

Hi emission is detected from Mrk 108 (see Figure [2]). The velocity field of Mrk 108 is 
distinct from its large neighbour NGC 2820 (Figure |2]^c)). The line profile of Mrk 108 is best 
fit with two components centered at 1406 and 1446 km s~^ with widths 52 and 117 km s~^. 
The component with centre velocity at 1446 km s^^ is contributed by NGC 2820. The 
emission from Mrk 108 is conti nuous to that from NGC 2820; hence we note that the mass 



from single dish observation by lThuan &: Martini (jl98ll ) is overestimated. The M(Hl )/Lk 



ratio is 0.41 and M{dyYi) / Lk is 2.86 adopting (V R/G) and the width of H I line is arou nd 



52 km s^^. A possible tidal dwarf galaxy is detected in the west (IKantharia et al.ll2005l ) of 
NGC 2820. Mrk 108 has experienced tidal interaction due to the large neighbouring spirals. 
Thus the situation is similar to Mrk 104 and Mrkl08 as a tidal dwarf galaxy cannot be ruled 



out. With masses being similar to those of tidal dwarf galaxy candidates (jWetzstein et al. 



20071 ) and being rotation dominated, the idea that a good number of BCDs in the groups 



might have originated in the tidal tails of larger galaxies is not an improbable scenario. 

The observed Hi emission of Mrk 1039 is about 1.1 times the size of the optical galaxy 
and the Hi disk is warped (see Figure [3t^a)). No warp is visible in the underlying stellar 
disk. The optical and kinematic centres seem to coincide. The velocity field shows a smooth 
gradient across the galaxy with some distortions in the western half. The higher resolution 
map in Figure EI^c) shows the higher column density regions on both edges of the galaxy disk. 
We note that radio continuum emission is coincident with the eastern half of the stellar disk 
hosting a strong H II region (see Figure |3]^f)) and a H I clump is also seen in the same region. 
High column density H I is also associated with the stellar component in the western half of 
the galaxy from which no radio continuum has been detected. The P-V diagram through 
the major axis (90° east of north; see Figure El^d)) also shows that the gas in the galaxy 
is mainly executing solid body rotation. The Hi line profile is double-horned (see Figure 
|8]). The fitted parameters are given in Table |3l The M(Hl )/Lk ratio for Mrk 1039 is 0.4 
Mq/Lq and M(dyn)/Li^ is 1.47. The width of the Hi line is 124 km s~^ 



The H I distribution and velocity field of Mrk 1069 are shown in Figure HI The H I disk 
(~ 185" diameter) is about six times the size of the optical disk (-D25 =~ 34") and higher 
column density gas is seen in the central and northern parts of the galaxy similar to the 
distribution of the radio continuum emission (see Figure St, d). Diffuse Hi is seen in the 
southern parts of the galaxy. The galaxy shows a smooth gradient in the velocity field with 
minor distortions. The optical and kinematic centre appear to be displaced. This has also 
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Fig. 8. — : The integrated Hi line profile of the galaxies Mrk 104, Mrk 108, Mrk 1039 and Mrk 1069 
are shown. The dots indicate the observed spectra and the solid line indicates the gaussian fits to the 
observed spectra. Note that while Mrk 104 and Mrk 108 show a profile consisting of more than one gaussian 
component, the spectra of Mrk 1039 and Mrk 1069 show the classic double-horned profile characteristic of a 
rotating disk viewed edge-on. 



been observed in other BCDs by iThuan et al.l (l2004l ) who attributed it to the energy input 
into the ISM through violent starbursts. The nearest galaxy, a dwarf, UGCA 052, is located 
around 16' (96 kpc) to the south western side of this galaxy but the two galaxies have not 
been classified in the literature as being associated. The M(Hl )/Lk ratio calculated for 
Mrk 1069 is 0.10 Mq/Lq and the width of line emission is 61 km s~^. The M(dyn)/Lx is 
about 1.99. 
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5. Conclusion 

The radio continuum observations at frequencies 610 MHz, 240 MHz, 325 MHz, 1.28 
GHz and 1.4 GHz with Hi observations and analyses of five blue compact dwarf galax- 
ies (BCDs) using the Giant Meterwave Radio Telescope (GMRT) are presented here. The 
Hi observations of four BCDs namely Mrk 104, Mrk 108, Mrk 1039 and Mrk 1069 have 
revealed their large Hi disks about 1.1-3 times the optical size of the galaxy. The typical 
H I masses range between 2 x 10^ — 10^ Mq. These values are typical of blue compact dwarf 
galaxies. We also detect a cloud close to the disk in Mrk 104 with no obvious optical coun- 
terpart and speculate that these might have influenced the current burst of star formation. 
Rotation is clearly noticed in the velocity maps of Mrk 108, Mrk 1039 and Mrk 1069. The 
Hi line profiles show multiple components. Mrk 104 has a Hi mass of 2 x 10^, rotation 
velocity is 69 km s"^ and M(dyn)/Lj^ ~ 1. Mrk 108 with total Hi mass is 1.6 x 10*^ Mq, 
M(dyn)/Li^ ~ 3, rotation velocity ~ 52 km s~^ and is situated very close to the galaxy NGC 
2820 (a spiral galaxy). Mrk 104 and Mrk 108 are similar with the nearest group member 
being a large spiral. A tidal origin cannot be ruled out for these two galaxies. 

All the galaxies are detected in radio continuum at 325 MHz and 610 MHz and the 
emission at this frequency is dominated by non-thermal emission. Out of the five galaxies, 
only two galaxies namely Mrk 104 and Mrk 1069 are detected at 240 MHz. 1 Zw 97 is 
detected for the first time in the radio; neither the NVSS survey nor the VLA FIRST survey 
detect it. Thermal and non-thermal separation is attempted for the four galaxies namely, 
Mrk 104, Mrk 108, Mrk 1039 and Mrk 1069. While the observed spectrum of Mrk 1069 can 
be explained by synchrotron spectrum absorbed by thermal gas at the lower frequencies, 
the observed spectrum of Mrk 104 could be explained as being due to combined thermal 
and non-thermal emissions. The observed spectrum of Mrk 108 is only due to pure non- 
thermal emission and the continuum emission in Mrk 1039 is explained as a combination 
of thermal and non-thermal emission with the non-thermal emission being absorbed at the 
lower frequencies by the thermal gas mixed in the region. We estimate a thermal fraction at 
1.4 GHz of ~ 80% and ~ 45% of the total emission for Mrk 104 and Mrk 1039 which is much 



higher than the ~ 10% generally seen in normal spiral galaxies (ICondonlll992l ). The emission 



in Mrk 1039, 1 Zw 97 and Mrk 1069 is confined to the H ll region and we note that the epoc h 



of star formation in Mrk 1039 is considered to be fairly young ~ 4 Myr (IHuang et al.lll999l ). 



The SFR estimated from the observed emission at 610 MHz results in values in the 
range 0.01-0.1 M© yr~^. These are si milar to the values of SFRs obtained for the individual 



star forming regions using Ha fluxes ( iRamya et al.ll2009bl ) in case of Mrk 1039 and 1 Zw 97 



whereas it matches with the SFR estimated from global Ha for Mrk 104. Emission detected 
at 610 MHz from Mrk 1039, Mrk 1069 and 1 Zw 97 is localized to the H ll regions seen in the 
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Ha images, which indicates that the SFRs obtained from the radio measurement represents 
the SFRs of a few individual H ll regions in these galaxies. On the other hand, the 610 MHz 
emission from Mrk 104 is seen to encompass the entire galaxy. 



Combining our results with those of iDeeg et al.l ( 1l993l ) — the only other study of BCDs 



in low frequency radio continuum — we find that on the average, the galaxies which are clas- 
sified as isolated galaxies tend to show a somewhat flatter observed spectrum as compared to 
galaxies which are classified as being in groups. This, we interpret as due to a larger fraction 
of thermal emission mixed with the non-thermal in the former case, as one expects for young 
localized starbursts. However a flatter injection spectrum of a young localized starburst 
region could also explain this. We further note that the starburst in a galaxy in tenuous 
environment is likely to be more localized. The flatter spectrum seen in isolated galaxies 
with localized star formation in the form of compact star forming regions and localized radio 
continuum emission, all suggest stochastic self-propagating star formation. 
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